Potential-flow models for channelled two-dimensional premixed flames 

around near-circular obstacles 
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The dynamics of two-dimensional thin premixed flames is addressed in the framework of math- 
ematical models where the flow field on either side of the front is piecewise incompressible and 
vorticity-free. Flames confined in channels with asymptotically-straight impenetrable walls are 
considered. Beside a few free propagations along straight channels, attention is focused on flames 
propagating against high-speed flows and positioned near a round central obstacle, or near two sym- 
metric bumps protruding inward. Combining conformal maps and Green's functions, a regularised 
generalisation of Frankel's integro-differential equation for the instantaneous front shape in each 
conflguration is derived, and solved numerically. This produces a variety of real looking phenom- 
ena: steady fronts (symmetric or not), noise-induced sub- wrinkles, flash-back events and breathing 
fronts in pulsating flows. Perspectives and open mathematical/physical problems are finally evoked. 



I. INTRODUCTION 

The propagation of flames through premixed gaseous 
reactants is governed by basic physical phenomena (sub- 
sonic and variable-density fluid mechanics, heat/mass 
transports, chemistry) that are now well understood in- 
dividually Ij. Also, modern computers allow their si- 
multaneous handling in a wealth of situations 0... if 
the flames are not of too large overall extent, L, or cou- 
pled with other phenomena (pressure waves, incoming 
eddies, radiant exchanges,...). This might suggest that 
the field is under control, but it's actually not quite so 
even for practical purposes, not to mention the theoreti- 
cal aspects. 

Direct numerical simulations must indeed solve the 
governing equations outside and inside wandering flames 
of thickness i <^ L, often a quite demanding task: The 
conduction-convection length £ = Dth/Uo = {DthtchY^^ 
based upon the fresh-gas heat diffusivity Dth, an over- 
all chemical time tch and the flame speed Uq of a flat 
(or nearly so) flame, usually is in the 10~^ m range for 
gaseous combustion at atmospheric ambient pressure pq 
whereas L typically is from a few 10~^ m to a fraction 
for laboratory or industrial burners; L/£ reaches several 
10^ in high-pressure combustion and hazardous explo- 
sions ■S'], rising to over 10^" for thermonuclear 'flames' 
in la-Supernovae [1, 0]. Up to 0(1) factors (Prandtl and 
Lewis numbers), £ also measures the range of counter- 
streamwise diffusions of momentum and species in gases. 
Furthermore, chemistry takes place in an even thinner 
internal layer of 0{£/Ze) thickness, where the Zel'dovich 
number of the overall reaction, Ze = 0(10 — 15), is an 
activation-to-maximum temperature ratio. 

A way out of the difficulties in many such situations 
is to idealise fiames as hydrodynamic discontinuities - 
fronts - separating chemically inert and inviscid fresh 
or burned media of markedly different densities, and 
equipped with a propagation law to fix the local burn- 



ing speed and with Hugoniot jumps, both amended by 
analytically derived corrections resulting from a finite £; 
the needed studies of local fiame structure are now essen- 
tially complete (see jfi-^]). Anyway, the number of nodes 
involved in numerical d-dimensional simulations of an in- 
terface embedded in the now chemistry-free flows still 
scales like {L/£Y, d =2 or 3, which often is the limiting 
step of simulations since L exceeds £ by far in most situ- 
ations. Restricting attention to the propagating front - 
the entity of primary interest - would decrease the num- 
ber of nodes down to {L/£)'^~^ and that of unknowns to 
1, but would require to theoretically derive an equation 
for the flame front itself, which has thus far not been fully 
completed. 

The difficulties that have been precluding a full theo- 
retical treatment of fiame front dynamics since the 1940s 
have three main origins: geometry (ies), nonlocal inter- 
actions and vorticity. Determining the moving front 
shape indeed constitutes a free-boundary problem, which 
is enough to bring about nonlinearities of geometrical 
origin even with linear field equations; the matter is 
further complicated by the often involved geometry of 
the domain where the front evolves (combustion cham- 
ber, obstacles, ...), especially in view of the difficulty 
to be evoked next that survives even when the combus- 
tion processes proper have a short range £ <^ L. Being 
markedly subsonic phenomena, premixed fiames cannot 
be described by a mere partial differential equation for 
their front, since low-Mach- number hydrodynamics com- 
bined with piecewise-uniform (but different) densities im- 
plies that all the fluid elements are coupled to the front 
motions instantaneously; the converse coupling cannot 
be neglected either, because a flame propagates relative 
to the fresh flow: the front evolution equation is also inte- 
gral in space, at least. As to vorticity, in general one can- 
not solve the fluid-mechanical equations for vortical flows 
even when piecewise incompressible and effectively invis- 
cid {L/£ ^ 1 implies a large Reynolds number), whereas 
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combustion invariably generates vorticity [9|; the task is 
not at all eased when an unsteady 'internal boundary' - 
the front - is wandering somewhere in the middle. 

To wit, the early works ('38 and '44) of Darrieus 
and Landau [ll| (DL) on the basic hydrodynamic flame 
instability, then on its 'ultra-violet' cut-off by local cur- 
vature effects had to consider nearly flat fronts, for 
only then could one solve the linearised burnt-gas hydro- 
dynamics and find the growth rate pertaining to a pre- 
scribed wave-number of wrinkling. The need to stay close 
to simple shapes also was felt when better accounting for 
the local physico-chemical flame properties [7|, or when 
studying the stability of near-spherical expanding fronts 
U2, f3]. It was also implicit in Sivashinsky's treatment 
[M] of saturation of the (then weak-) DL instability by 
weak nonlinearities in the limit of small density changes, 
and in extensions thereof all addressed nearly-flat, 
or -cylindrical, -spherical, fronts. 

Until recently the sole exception to the 'simple-shape 
& weakly-nonlinear' constraint, allowing for large front 
wrinkles while retaining nonlocal aspects of front-to- front 
interactions, was the equation proposed by Frankel [l^ 
as a potential-flow extrapolation of Sivashinsky's to arbi- 
trary amplitudes of wrinkling. The latter equation, orig- 
inally obtained ,14] from a flow that happens to be po- 
tential to leading order for small density changes, turned 
out to also govern the shape of steady flat-on-average 
wrinkled flames when the next order (hence some vor- 
ticity effects) is retained [13, lisl: retaining up to three 
more orders for steady fronts leads to a gener- 

alisation known as the Zhdanov- Trubnikov equation [23| 
with similar properties..., just like potential flow models 
do [llj. As shown later on here, integral formulations 
via Frankel-like equations show good prospects to handle 
the complicated geometries of the flame front along with 
the influence of the boundaries and the nonlocal inter- 
actions; it may ultimately prove to constitute a mostly 
convenient way of handling thin flames even wheii vor- 
ticity is accounted for: as a perturbation like in 22;] or, 
better, on properly adapting the non-perturbative meth- 
ods of [1^,11^ (yet to be done). 

Recent advances indeed made it possible to relax the 
restriction to weakly nonlinear shapes while still retain- 
ing vorticity (0, Hj], and references therein), in simple 
large-scale configurations: two-dimensional spatially 2L- 
periodic flow- field and front, flat-on-transverse- average 
flame, smooth spontaneous growth of disturbances since 
the remote past, piecewise-incompressible fluids. It was 
shown that the task of solving the Euler equations can 
then be bypassed when only the flame-front evolution is 
sought, which resulted in closed equations for the front 
shape and potential velocity components defined along it 
(as opposed to 'in the bulk') whatever the front wrinkle 
amplitudes and the expansion ratio. Such equations are 
necessarily nonlocal in space (to account for the nonlo- 
cal potential interactions like in Frankel's) and in time 
(to keep track of the flame history, stored 'space-wise' 
as a vorticity field), except when the burnt gas flow 



is quasi-steady; also, these equations are not yet fully 
ripe for numerical treatments, especially for unsteady 
fronts, and are still restricted to such simple bound- 
aries as straight channels. Except in rare 'no-boundary' 
situations (e.g., unconfincd expanding flames), labora- 
tory measurements or specific appliances do involve more 
complicated boundaries (and conditions along them), 
whose part precisely is to affect the flame behaviours in 
a (hopefully-) controlled way, e.g., to keep them where 
required. How to incorporate such boundaries remained 
to be done, which a priori is a rather inconvenient con- 
straint for comparisons with theories, since flames of in- 
finite lateral extent are not easy to study experimentally. 

Recall that the equation derived in [l^ for steady 
fiames in straight channels (or of infinite lateral extent) 
takes the form of an integral equation for some poten- 
tial velocity field with vorticity-affected jumps across the 
front; these resulted from a delicate simplification of the 
solution (stream-function) to a Poisson equation, so as 
to isolate its vorticity-containing downstream contribu- 
tions near the front [1^, 113 . The nonlocal equation de- 
rived by Frankel for d = 2 potential-fiow models also 
results from a Poisson problem for a ( complex- valued) 
stream function, like the nonlocal equations in [isl . [23l | 
yet only with concentrated non-negative sources along 
the front, whereas actual flows also have a distributed 
'source', vorticity, of zero transverse average. As such 
it has a family likeness with Kazakov's equation (i.e., 
in steady cases) [isj and shares important ingredients 
with it; most notably the Green's function of Laplace's 
equation that is long-ranged in 2d and hence constitutes 
a convenient vehicle to encode geometrical properties of 
boundaries: shape of channel walls and obstacles to keep 
fronts where needed, for example. 

Obtaining Frankel-type equations for fronts in the 
presence of obstacles, taking up the numerical treat- 
ment/difficulties, and displaying the associated front 
properties may thus prove of interest per se, as a flame- 
like dynamical system able to cope with more realistic 
geometries (e.g., multiply-connected fluid domains) than 
in [i^, while accounting for nonlocal self-interactions, 
and as a prerequisite for vorticity-containing formula- 
tions; comparisons between the two will also help un- 
derstand and quantify the specific influences of vorticity 
effects: insofar as spontaneous front evolutions (or occur- 
ring on the same time scale if forced) are concerned, the 
influence of vorticity indeed turned out to be only quan- 
titative (weaker-than-potential hydrodynamic instability, 
different constants in evolution equations. . . ). Potential- 
fiow flame models with nontrivial geometry of boundaries 
are what the present contribution is about, in a few d — 2 
configurations accessible to complex-variable methods, 
conformal mappings in particular: the fiexibility of such 
tools gives access to interesting domain shapes without 
affecting the Neumann conditions at impenetrable walls 
or the Laplace equation in the bulk; also. Green's func- 
tions trivially transform. Nonlocal front self-interactions 
also are accounted for in a natural way. 
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Our (admittedly restricted) present scope thus only 
concerns two aspects -geometry (ies) and space- wise non- 
local interactions - of the complete, as yet unsolved, prob- 
lem. We are fully aware that omitting vorticity can possi- 
bly miss some aspects of flame dynamics, and only qual- 
itative agreement with actual fronts is expected in gen- 
eral, rendering quantitative comparisons with full-blown 
direct numerical simulations and experiments still prema- 
ture; there could be more than that on occasion, however, 
since vorticity and its conceivable consequences only lit- 
tle show up in some important instances (see [IJ]) and 
are then likely amenable to perturbative treatments like 
in [23. Moreover, given the structure of the vorticity- 
affected equation derived in [l^ , the tools and results of 
the present contribution will also be of interest to tackle 
the full problem. 




negligible throughout. Hence chemical, or transport- 
related, front /boundary interactions like wall ignition or 
anchoring [23 are not addressed here. The flat-flame 
speed C/q is assumed known from a separate analysis (e.g., 
see [ij), as is the Markstein length C introduced in equa- 
tion ^ below. Density is taken as unity over the fresh 
gas side and 1 / E < 1 in the burnt medium, whereby the 
velocity tt, with components (u, v) in Cartesian coordi- 
nates X = {x,y), Fig.[T] is solenoidal on either side: 

V-M = 0. (1) 

Outside of the front, the no- vorticity assumption inherent 
to Frankel types of models implies u = Vcp for some 
potential (/?, whence ip itself is harmonic: 

VV = 0. (2) 

Even though potential flows can accomodate viscous ef- 
fects at interfaces [2^, the simplest jump relations com- 
patible with a nonzero mass flux across the flame [l[ are 
adopted here: the tangential gradient n • r of and the 
potential itself are considered continuous at the front, 
whereas the normal velocity satisfies 

[m • n] = u ■ n|+ — u ■ = {E — 1) U() a-{x, t). (3) 

on account of mass conservation and density jump [l6| . 
Here n is the unit normal to the front pointing to the 
burnt side, and t is time. The right-hand side of ^ 
measures the rate of 'volume production' (per unit front 
area) associated with the density drop from 1 to < 
1. Finally, the boundary conditions associated with ([2]) 
include 



FIG. 1: Sketch of a d = 2 flame about a central obstacle, and 
coordinate system. 

The paper is organised as follows. The regularised Po- 
tential Models and their general integral form are intro- 
duced in Sec. 2. The selected method of numerical in- 
tegration is presented in Sec. 3, in the context of free 
flames confined by parallel channel walls. Section 4 
addresses channelled flames settled counter-streamwise 
about a round central obstacle (Fig.[T]) in high-speed 
flows, the related case of front positioned near lateral 
bumps being examined next, in Sec. 5. Flames fed by os- 
cillating fresh flows are taken up in Sec. 6. The paper 
ends up with conclusions, perspectives and open prob- 
lems (Sec. 7). A short Appendix analyses slender flames. 

II. POTENTIAL-FLOW MODELS 
A. General formulation 

The flow is henceforth assumed two-dimensional, 
piecewise incompressible on each side of an infinitely thin 
flame- front. The internal time ^/C/q ~ ich is considered 



u- rij =0 (4) 

along the j-th impermeable boundary with local nor- 
mal rij. In channel geometries like in Fig.[Tl u also 
is prescribed to be uniform and parallel to the y-axis, 
(m, v) = (0, V) ai y — —00 (incoming fresh-gas flow); for 
free propagation and steady flames (if any) V = Vp is 
a kind of eigenvalue to be found as part of the solution, 
whereas for flames near obstacles like in Fig.[Tl V = Vi„ is 
a prescribed, high-enough (I^n > Vp) injection velocity. 

Given a 'charge-density' Uo(j{x, t), ([I])-(|4|) are sufficient 
to compute u{x, t) pertaining to a presumed flame shape 
and, in particular, the fresh-gas velocity u-n\^ just ahead 
of its front. To determine the latter's evolution, its nor- 
mal velocity T) is prescribed to satisfy: 

n • u|_ - X> = [/ocr(a3,i) > 0, (5) 

where the (scaled-) burning velocity a may depend ex- 
plicitly on X. This could help one model 0(1) variations 
in JJoj without noticeable change in the densities, caused 
by 0{\/ Ze) <^ 1 fractional modifications in mixture com- 
position; or in models where a{x) would account for an 
effective local burning speed, determined by averaging 
over small(er)-scale wrinkles. Like in Frankel's original 
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work [TgI we shall assume a{x,t) = 1 throughout the pa- 
per, even though the cr(-, t) function will be kept in most 
formulae, for completeness. To simplify the subsequent 
formula we shall henceforth adopt units where Uq = 1. 

With (7 = 1 in ([J])©, the problem is dynamically ill- 
posed, however, as is seen on considering a steady flat 
front {Vf — 1) located at y—Y— const, and perturbed 
by infinitesimal wrinkles proportional to exp{ikx + wt); 
routine Fourier analysis of the linearised ©-© yields 
{E — l)|fc|/2 for the wrinkle growth rate w. Because 
the above formulation is of first order in time and au- 
tonomous (one may then everywhere omit t as argument 
to save space), it has no inertia/memory effect and tu/j/cl 
exceeds what the DL analyses gave, implying a stronger 
mechanism of spontaneous wrinkling: the DL growth rate 
w > satisfies 2m \k\ = P{E - 1) - {1 + l/E)w'^. As 
expected on dimensional grounds (only one way to get 
a left o right invariant growth rate w from a flat-flame 
velocity Uo{= 1), a wave number fc, and two densities), 
both predict 137 ~ with n?/ |/c| > depending on E 
only and vanishing at E = 1. At any rate, w — ?> +00 
as |fc| — T' -l-oo when E > 1, thus rendering the dynamics 
ill-posed, and an expedient way to regularise the prob- 
lem is to introduce a Markstein cut-off length £ > @ , 
replacing ^ by 

u-n\^-V = a{x,t)-CC, (6) 

where C is the local front curvature. C is proportional 
to the actual conduction/convection flat-flame thickness, 
£, but exceeds it by a large enough numerical factor (15, 
say) that one may retain curvature effects in a thin front 
model. The new w{k) deduced from ^ reads zu = {E ~ 
l)|fc|(l — \k\/kn)/2 when (T{x,t) = 1, with the neutral 
wave number: 

kn = iE-l)/2C. (7) 

Besides restoring well-posedness, a related virtue of the 
ultra-violet cut-off is to smooth the otherwise sharp front 
cusps (crests) that would form. We stress that in the 
present thin-flame formulation, as well as in those re- 
taining vorticity (isl . [23j , all the physico-chemical prop- 
erties of the mixture are lumped in the fresh-to-burnt 
density ratio E > 1, in the flat-flame velocity Uq (taken 
here as unit of velocities), in the fresh gas density (se- 
lected as unit of densities) and in the Markstein length 
(more generally, in the local propagation law). To ease 
comparisons with experiments that would have access to 
kn, the Markstein length in ([6]) could be expressed as 
{E-l)/2kn. 

B. Integral form 

It is convenient to introduce the usual notation: w = 
u + iv for the complex velocity field, $ = ip + iijj for 
the complex potential, ip being the stream-function con- 
jugate to If, z = X + iy for the complex variable in 



physical plane; and over-bars for complex conjugation, 
so that w = u — iv = d^/dz. We also define the com- 
plex forms of X = (TxTTy) and n as r(s) — + iTy and 
n(s) = 'iT{s), respectively. These are here evaluated at 
a point labelled by the arc-length coordinate s along the 
front z{s) = x(s) + iy{s), oriented in such way that the 
burnt gases lie on its left as s increases. 

Let G{z, z') be the complex potential associated with 
a unit source of fluid located at z' in the fluid domain, 
satisfying 2'kG{z, z') — >■ ln(z — z') for z — >■ z' and the slip 
condition Im(G') = i^j{z') along the jth impenetrable 
boundary. An important property of the Green's func- 
tion G(z, z') in such channel-like domains as in Fig.[TJ 
without or with (a finite number of) impenetrable obsta- 
cles between the lateral walls, is worth noticing. Consider 
a streamline along which Im(G') = ipLiz'), issued from z' 
and ultimately reaching the leftmost channel wall where 
it splits to follow the wall to ±?oo; it is readily seen to 
be unique, as is the other one Im(G) = V'fi(^') ultimately 
meeting the right wall, and doing the same there. The 
mass fiow rates to z = ±ioo are equal to \tPr — tPlI and 
hence both equal 1/2 since they sum up to the unit yield 
of the source. This property holds whatever z', and is in- 
variant by conformal mappings from the physical plane 
to another, since the mass flux crossing a curve is [27l |. 
Thanks to the linearity of ^ one may superimpose vol- 
ume sources along the front and ensure that ^ is ful- 
filled. The resulting cumulative complex velocity reads 

wJsources = {E - l)^ j G{z,z') a{s') ds' , (8) 

where the integral is over the whole flame length, z' = 
z(s'), and f7(s) is shorthand for a{z{s)). As emphasized 
above, exactly half the total yield (E—l) J cr{s')ds' of the 
sources along the front flows to z = ±ioo. At z = —zoo 
this produces a uniform flow. To compensate for it and 
satisfy the condition v = Vln at ?/ = — c», a supplemen- 
tary flow Wsup(-z) must be added to Wsourccs(-z), such that 

Wsup(-joo) = -i / o-(s') ds' - iVin- (9) 

This Wsupiz) is analytic across the front, because (|3]) 
is already accounted for by ([5]), and must also satisfy 
the slip condition ^ along all impenetrable bound- 
aries. Wsup{z) /wsupi—ioo) will only depend on the ge- 
ometry of the channel and/or obstacle(s) under consid- 
eration and is closely related to the conformal map from 
the physical channel to an auxiliary straight one (en- 
dowed with cuts for each obstacle(s) present); in partic- 
ular, Wsupiz) = Wsupi~ioo) for straight channels. With 
w = Wsourccs+Wsup in principle available in the bulk of the 
fluid, one may specialise it to the entrance of the front, 
then substitute in ([5]) to obtain an evolution equation for 
the front z{s,t), e.g., see equation ^5]) below. 

Before proceeding any further we note a general 
(necessary-) test of steadiness, obtained on integrating 
Eq. Il]) over the whole fresh gas domain; by ^ (and 
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iC = d\n{T{s))/ds, J C{s)ds = here since t(s) is the 
same at the right and the left front ends), this gives 



(10) 



where Y{t) is the x- averaged front location along the y- 
axis: Steady flames have 14n = / cr^s) ds/2L, whereby 
cr(s) = 1 implies that the speed of a steady flame as a 
whole relative to the upstream fresh medium simply is (in 
units of the flat-flame speed) the front-length to channel- 
width ratio; for unsteady flame patterns, Vpit) may in 
fact be defined as J a{s,t) ds/2L. Only relying on the 
local propagation law © and overall mass conservation 
in the fresh medium, (jlOp would still hold true even with 
vorticity retained. 



III. FREE FLAMES IN STRAIGHT CHANNELS 
A. Green's function and evolution equation 

To the best of our knowledge, propagations of 
markedly wrinkled front have so far been studied nu- 
merically by means of potential flow models only for ex- 
panding fronts [^,[13, i.e., the particular situation that 
Frankel's equation envisaged, and for Bunsen types of un- 
confined flames attached at points [11] (or along a circle 
j^ol], for d = 3) on an injection line (or plane); we consider 
here the case of free propagations along straight channels 
with parallel walls a.t x = ±L without any obstacle in 
between, which had so far not been investigated. The 
Green's function pertaining to Neumann conditions at 
the walls, then denoted Go{z, z'), is the familiar [27l[30|: 

2nGoiz,z') ='ln(sin-^(z-z') cos-^(z + z')) , (11) 

where the symbol =' means 'equal up to additive func- 
tions of the primed quantities only', here z'; these are 
indeed annihilated by dz, and hence would play no role 
in ([5]). The cosine term in (ITU) accounts for the images 
in the lateral walls, and Gq is 4L-periodic in z. 

Carrying the differentiation needed in ([S]) and special- 
ising the result to the entrance of the front yields the 
contribution to fresh gas velocity ahead of it induced by 
the sources: 

WsourcosI- ^ l^(T{s)+j-Wo{z{s),z{s'))^^^ds' 

(12) 

the first term in the right-hand side follows from the 
Plemelj-Sokhotsky formulas for the limiting values of in- 
tegrals of the Cauchy type [27] as z approaches the inte- 
gration contour (the front, here). In p^ . and subsequent 
formulas, {Tr/4:L)Wo{z, z') will represent 2TrdzG{z, z'); 
explicitly: 

Woiz, z') = cot -^{z- z') - tan -^(z + Y). (13) 



As for the supplementary flow Wsup needed to satisfy the 
conditions at y = — oo, it is uniform in the present con- 
figuration, hence is given by (|9]). We now have all the 
ingredients to write 



n • M_ = Re {— i n(s)} ^Vf + ^ ^ g(s) ds 
-|-^^^Re| I n{s)WQ{z,z')a{s')ds' - \La{s) 



(14) 



To obtain (jl4D we used the fact that the scalar product 
of any two vectors a — (a^, ay) and b — (b^, by) is a b — 
Re(a6), with a = {qx -\-iciy) and similarly for b. Whereas 
It • t]_ requires principal parts, the above integrals are 
ordinary ones: since tan(a) = a + a'^/3 + . . . for \a\ — >■ 
0, Re[n(s)cot(7r(z - z')/4:L)] -2LC{s)/tt as s' ^ s 
and the integrand in ([T4|) is finite if z{s) does not lie 
on one of the boundaries; if it does, Ke[n(s)Wo{z, z')] — 
-ALC{s)/t: at s' = s. 

The evolution equation for freely propagating flames 
in straight channels can now be written as: 

ReJn(s)(TZ;_ - 9(z(s, t))] = 1 - £C(s, t), (15) 

with Re(n(s)i(7_) given by (HH). 

For infinitesimally wrinkled front, z{s) = x{s) + 
iF{s,t), dxF — > 0, one may write t(s) — 1, x{s) — s, 
Vp = 1 = Vin, C — dxxF, and linearise the integral term 
of (fH]) after replacing it by a principal value because 
the linearization procedure is not uniformly valid in z 
and z' , and extend the integrand to a 4L-periodic func- 
tion. Integration by parts then gives dtF ~ CdxxF + 
\{E— l)H{—dxF}, where the 4L-periodic Hilbert trans- 
form H{-} has zA;exp(ifca:)} = \k\ exp^ikx). For the 
wavenumbers allowed by ([4]), the same dispersion rela- 
tion 'uj{k) = {E — 1) \k\ (1 — \k\ /kn)/2 as before is thus 
recovered. 

If now < - 1 < 1, ru < 1 for ah \k\ /fc„ < 0(1) 
and a weak nonlinearity can counteract the instability 
[l^ . With knL — 0{1) assumed, geometrical effects 
merely add the presumed-small —{dxF)'^/2 ^ {knF)'^ 
to the above dtF ^ voF ^ [E — l)knF <^ 1, whereby 
dxF ^ knF — 0{E— 1) <C 1 is small, as anticipated; this 
nonlinearity stems from the near-unity reciprocal cosine 
(1 + {dxF)^)^/^ = 1 + {dxF)^/2 ... of the angle be- 
tween n and the mean direction of propagation (y-axis, 
here) that enters the definition of the normal front ve- 
locity V = dtF/{l + (9^F)2)i/2, Ultimately, the leading 
order equation for F(x,t) reads: 



dtF+1- 



(dxF)' 



E-l 



H{-dxF} + V,^, (16) 



which is Sivashinsky's [i4|: a 'synthesis' of the lin- 
ear result and of an eikonal equation driven by the 
geometry- induced {dxF)^/2. To second order in the 
{E — 1) expansion, the equation for F in steady cases 
keeps the same form, up to (removable) .E-dependent 
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coefficients [T3, [II], and the left-hand side only ac- 
quires a single additional nonlinearity proportional 
to —JH{—dxF})'^ when two more orders are retained 
[isl [igj . e.g., stemming from the —{dxF)'^/8 term of 
the expanded (1 -I- {dxF)^y^^; the equation for F{x) 
then becomes of the Zhdanov- Trubnikov (ZT) type [20l |. 
At any rate, equations with both nonlinearities can be 
solved in terms of the same elementary functions (and 
coupled nonlinear ordinary differential equations) via 
a pole decomposition technique [3l|, [s^. Some issues 
relevant for the present work can be summarized as 
follows: 

(i) Thanks to the stabilizing nonlinear-term of geomet- 
rical origin, generically leads to steady bi-coalesced 
fronts (sharp maxima of F where dx = B,t each wall, and 
only there) when integrated with the Neumann boundary 
conditions implied by steadiness of course requires 
that Via = Vf = 1 + / {dxFY dx/AL in accordance with 
the small- front-slope form of (ITUl) . (ii) Several stable bi- 
coalesced solutions with the same Vp, 0{4,L{E —1) / AirC) 
in number, exist for the same mixture properties and 
channel width 2i [H Hi], (iii) For channel widths 2L 
that noticeably exceed the half neutral wavelength A„/2, 
A„ = 27r/fc„ = A'kC/{E — 1), the stable steady solution 
profiles comprise inner regions of 0(A„) ^ L width ad- 
jacent to the walls where the curvature term and C must 
be retained to ensure dxF{x — ±L) — 0, separated by an 
outer zone where one may omit CdxxF; the outer profiles 
are the restriction to —L < x < +L of 4L-periodic pat- 
terns that have sharp cusps and satisfy the weaker bound- 
ary conditions dxF{zLL — x) = —dxF{zLL + x), and are 
available analytically [3l|, [SJ] . (iv) Many unstable steady 
solutions exist, and their number grows with 4L/A„ much 
faster than linearly [sl, HI], (v) The slightest addi- 
tive random noise in the right-hand side of (|16p makes 
the flame execute random transitions between metastable 
states with small 0{{E — 1)^A„/L) fractional variations 
in front length (33| and instantaneous a;-averaged speed, 
which can be partly understood in terms of extra sin- 
gularities (poles) of dxF sporadically implanted by the 
external forcing [3^ . (vi) These properties hold for both 
(IT51) and its ZT generalisation, i.e., persist rather deep in 
the {E — 1) expansion. 



B. Numerical free flames 

Essentially because (1 + (dxF)'^)^/'^ <l + {dxFf/2, it 
was not obvious from (jl4p that an unexpanded stabilising 
geometry factor alone would be able to saturate hydro- 
dynamic instability at larger scales than A„ when E—l = 
0(1), though several qualitative arguments suggested it 
could. Firstly, this stabilising influence has a simple in- 
terpretation: if the wrinkle-induced velocity component 
in the local mean direction of propagation gets too high, 
the front can avoid being blown-off by rotating locally 
to ensure that ^ be satisfied (like flames above Bunsen 



burners, or the stabilized fronts to be considered next), 
which brings about the cosine {l + {dxF)'^)~^/^ and gener- 
ates curvature-smoothed crests whenever adjacent front 
pieces rotate in opposite directions; this will likely survive 
for large amplitudes of wrinkling. Secondly, an equation 
of the form 

(1 + {dxF)'y/^ = ^_lH{-dxF} + Vf, (17) 

where the reciprocal cosine (1 + ((?j;F)^)^/^ is not ex- 
panded, can be recast (by squaring it) into the A„/L = 
0+ form of a ZT equation that, like the £/L{E-l) = 0+ 
steady outer version of with T^n = Vp, possesses 
bounded solutions F{x)y^ with sharp cusps at the walls 
(but a bounded Vp) [33]. Finally, expression giving 
n • m]_ — Vp for a presumed F{x) is 'less severe' than the 
driving term (E — l)H{—dxF}/2 of (fTB|) . in the sense 
that the contributions of the integral over s' in (fT4)) are 
comparatively less sensitive to large amplitudes of wrin- 
kling; the reason is that the trigonometric functions in 
P3p are bounded when their arguments get large imagi- 
nary whereas the linear operator H{-) with a real cotan- 
gent as kernel does not 'saturate' for \y{s') — y{s) \ » L, 
which likely results in a milder source of wrinkling than 
in (dll). 

To determine the front wrinkle evolution/shape(s) 
when E — 1 = 0{1) and significant nonlinear effects set 
in, we resorted to a numerical integration of (ITil)(|15l) . 
To this end, the integral featured in (fH]) was discretised 
at nodes sq = 0, si, . . . , sat, and computed by the trape- 
zoidal method, account being taken that the integrand 
in is continuous at s' = s (see above Eq. p^ ). The 
local front curvature C(s,„,t) is calculated by passing a 
circle through z{sm-i,t), z{sm,i), z{sm+i,t), which also 
gives n{srmt) then T{sm,t). The normal node displace- 
ments D{sm,t)n{sm,t)St resulting from (ITHlp^ enable 
one to locate the flame at t + St. To ensure sufficient 
resolution while still avoiding the numerical stability cri- 
terion 2£St < (<5s)^ from being violated in this explicit 
scheme, the arc-length step Ss was kept within bounds 
[<^Smin, ^Smax] choscn to compromisc between CPU cost 
and accuracy; typically Ssmin — L/200, (5smax = 4(5sinin, 
N = 300. Monitoring the s-grid requires some re- 
meshing/interpolation, which brings about small-scale 
numerical noise. As checks, we verifled that stable flames 
resulted when the parameter 

i, = 7rC/{{E-l)L), (18) 

exceeds unity; l/v is a Peclet number based upon the 
neutral wavelength (A„ = 27r/fc„) and that of the wrin- 
kle {— 4L, in absence of a; —x symmetry). We also 
checked that < i? — 1 <C 1, fc„L = 0(1) give flame 
shapes and speeds that are compatible with Sivashinsky's 
equation. Before proceeding to results, we note that free 
flames enjoy a Galilean and translation invariance along 
the y-axis. Two main ways thus exist to compute the 
speed of steady patterns: adjust V^n until the x-average 
front location, Eq. (ITUl) . also stays steady; or set T^n = 1 
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FIG. 2; Steady free fronts numerically obtained for E — 5, 
L = 2, £ = 0.25; see the main text, Sec. IIIB. 

and wait until the front recedes towards y = — oo at the 
constant speed 1 — Vp < 0. The second, easier, method 
was adopted. 

A first noteworthy numerical fact is that steady free 
flame patterns do exist when E — 1 = 0(1), as is evi- 
denced in Fig.EJfor E = 5, L = 2, C ^ 0.25: due to the 
small £/L ratio involved in this case, curvature effects 
are mainly felt (Fig. [2]) at the front crests, which is not 
true any longer when the parameter v defined in (jlSp 
approaches unity from below. 

Like with Sivashinsky's equation [H, [s^l when Neu- 
mann conditions are imposed at the lateral walls, the 
flame shape generically has no a; — >■ —x symmetry. Yet 
symmetric solutions also exist if the channel is wide 
enough. Figure [2] settles the point, the symmetry of 
the solution labelled 1.47 being then imposed by numer- 
ically solving l|14p (fT5|) over a half-channel with the same 
C, then completing the front-shape function according to 
F{—x) — F{x); otherwise, the central extra crest in Fig. [2] 
would ultimately join x — ±L whenever shifted to the left 
or right by the symmetry-breaking imperfections (numer- 
ical, or via initial data). The above mentioned procedure 
to generate unstable multi-crested fronts for wide chan- 
nels can be extended by concatenating replicas of halved 
periodic patterns of a smaller minimal wavelength A than 
2L (with 4L/A = m = integer), but such flames with sev- 
eral crests would also be unstable against slight lateral 
shifts or compression of any of them; the unstable flame 
fronts labelled 1.47 (to = 2) and 1.30 (to = 3) in Fig.H] 
were obtained in this manner. The solutions labelled 1.49 
and 1.58 were both obtained by the time marching nu- 
merical method without the 'replica-trick', however, and 



hence represent two different (meta-)stable solutions of 
the same problem in the same operating conditions! 

Thus, the picture is somehow reminiscent of that de- 
duced for (S- 1) < 1 from HH) and its ZT higher-order 
generalisation, when endowed with Neumann conditions 
(see end of the preceding subsection): many unstable 
steady states (at least 0(4L/A„) in number, thanks to 
the 'replica trick'), several meta-stable steady fronts; 
to which one may add the sensitivity to noise when 
4L 3> A„. However, no such steady shape as that labelled 
1.49, requiring well chosen initial conditions, could so far 
be obtained by the pole-decomposition method, and it 
possibly is of a new type: contrary to the ones belong- 
ing to Vf = 1.30 or 1.47 in Fig.[2J it has no apparent 
symmetry. 

The solutions in Fig. [5] have different front lengths and 
effective speeds, Vp = 1.58 or 1.49 for the asymmetric 
patterns, and Vp = 1.47 or 1.30 for the A = 2L or 3A/2 = 
2L symmetric ones (and still smaller Vp — 1 > for 
771 = 4, 5, . . . as long as AL/m > A„), max(yF) exceeding 
what direct numerical simulations predict, Vp ~ 1.20 for 
E = 5 (see [l^ and the Refs. therein); this likely is 
related to the present larger linear growth rate w{k) at 
|fc| <C fc„ evoked above (O. Note that quite different 
flame patterns may have nearby V^s. Accumulated runs 
revealed that Vp — 1, ^Inax — ^min and CC{x = L) increase 
with E\ this was expected because all vanish at, and grow 
in a neighbourhood of, = 1 just like the linear growth 
rate vj{k <C fc„)/ \k\. 

IV. CHANNELLED FLAMES SETTLED NEAR 
OBSTACLES 

A. Conformal map, Green's function 

As mentioned earlier, the first step to generalise the 
evolution equation ((Ti)) to the situation sketched in Fig.[T] 
is to compute the Green's function G{z,z') belonging to 
it. The easiest way to it is to tailor a conformal mapping 
Z{z) from the fluid domain in physical z-plane (between 
the lateral walls and exterior to the centred obstacle) to 
some Z-plane where G, then noted G{Z,Z'), is easily 
accessible, viz.: 

G{z,z')^G{Z = Z{z),Z' ^Z(z')). (19) 

For any single obstacle the Riemann mapping theorem 
[27I I guarantees the existence of such a map from the 
doubly-connected fluid domain to the annulus q < \Z\ < 
1, for some < g < 1 (Fig. [3]). Ideally, one would 
start from a circular obstacle \z\ < r of given radius 
r < L, then determine Z{z). We proceeded the other 
way around, first noting that the map 

TT Z 

Z(z)=tan— (20) 

transforms the two walls x — ±L into the circle \Z\ = 1, 
X = +L (or —L) corresponding to Re[Z] > (or < 0), 
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FIG. 3: Annulus g < |^| < 1 in auxiliary Z-plane. The 
streamlines on the left are for a cold flow injected at Z — —i 
(potential '^{Z, q) defined in (|29|) 'l and removed a.t Z = +i. 

and z — ±ioo to Z = ±i; next, that a small enough 
domain of 0{qL) size around z = is mapped by (^0)) 
to a domain of nearly identical shape if q ^ f, as a 
result of tan(a) = a + a-^/S + . . . for \a\ — > 0. We thus 
postulated that the obstacle image by (l20l) exactly is the 
circle \Z\ = q < I, its shape in physical domain being 
obtained through (I^Dl) as 

z = - arctan((je''^), -tt < T < tt, < g < 1. (21) 

q 

Since its aspect ratio is ry/r^ — atanh(g)/ arctan((j') = 
I+2g^/3 + . . ., Fig.[T], a fairly-circular physical obstacle is 
obtained whenever q < 1/5, i.e., r = ALq/ir < L/3.93 . . ., 
which is amply sufficient for our present purposes. 

The analytically invertible map d^H]) to the annulus 
q < \Z\ < 1 also is such that G{Z,Z') is accessible in 
closed form: 

2.Giz,z')^' ,jnziz',q)nzw^q)\ 

^ ' \ P{iZ,q)P{-iZ,q) ) ^ ^ 

in terms of a 'loxodromic' function P{h,q) defined [ssl ] 
as: 

P{h,q) = {l-h)\l{h,q), (23) 

OO 

U{h,q) ^ \{{l~q^^h){\~q^^/h). (24) 

Using (23]) (HH), it can be shown that 

P{q^h,q) = Pil/h,q)^-P{h,q)/h (25) 

P(M) = P(h,q) (26) 

These properties and Gauss's triple-product identity 
[38l ] help one show that the above G{Z, Z') has a Z- 
independent imaginary part along \Z\ = q and along 



\Z\ = 1 except for jumps at Z = ±i where two sinks 
reside. These have unit weights, which seemingly contra- 
dicts the remark made above ([8]): the north/south unit 
sinks sit on the boundary, however, whereby each con- 
tributes only —1/2 to the overall mass balance in the 
annulus (Fig. [3]). 

Because needs be differentiated in expression 
of {u — TO)sourcos, wc introduce the notation: 

IC{h,q) = h^lnU{h,q) 

I °° 1 

" ^h^^^^^^q-^^ + q^'' -h~l/h ^^^^ 

This IC{h,q) is to be ultimately evaluated at h = Z/Z', 
h = Z'Z' and Z = ±i, and d^U-dSZ]) clearly displays 
the infinite series of images oi Z = Z' in the two circles 
\Z\ = 1 and \Z\ = q {at Z = Z'q^'' and Z = -g^V^, 
k = ±1, ±2 . . .) that are all needed (together with those 
of Z = zLi) to fulfil the Neumann conditions along the lat- 
eral walls and the obstacle, and mass conservation; since 
the inverse z{Z) of (l20l) is defined up to integer multiples 
of 4L, each point in the Z-plane corresponds to an infi- 
nite series of pre-images in the physical plane, 'beyond' 
the lateral walls. Note that /C(/i, 0) = = /C(±l, q). 



B. Supplementary flow 

The second contribution to w = ^sources + wJsup also 
is easily expressed in terms of the P{h,q) function as 
d^sup/dz, with 

$sup = 4L (^—^ f a{s') ds' + Mn) ^{Z, q), (28) 



^{Z, <?) = ^ In (P{-Z/i, q)/P{Z/i, q)) . (29) 

The potential 4'(Z, g) indeed corresponds to a pair of 
weight-two source {Z — ~i) and sink {Z = +i) sitting on 
the circle \Z\ = 1; some associated streamlines in Z-plane 
are displayed in Fig.|3] and correspond to an inert gas 
flowing past the obstacle. On combining P(/i, 0) = (1 — ^) 
and (|20[) . one can verify that the above supplementary 
flow resumes the uniform one found in subsection 3.1 
when g — > (at fixed z ^ 0), as it should. 



C. Evolution equation 

The evolution equation has the same form as (IT5|) . the 
only difference being that the fresh-gas fiow velocity nor- 
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mal to the front now is given by: 

n ■ M_ n • Uo_ + Ke[n{s){Z H )JC{iZ, q)] 

2 Z 

" (^ + ^ ycT(s)d.)ReK.)(Z + l)/CHZ,g)] 

I ^(^')^ReKs)(Z + i)/C(|,<z)] 
+ ^ / a(.')^ReK.)(Z + l)/C(ZZ',g)], 

(30) 

where n-Uo_ formaUy coincides with the right-hand side 
of (ITil) . yet evaluated in terms of Via (instead of Vp) and 
the current unknown z{s,t). Recall that )C{-,q) vanishes 
at g = (no obstacle), as does Z +1/Z = 2/ sin(7rz/2L) 
when the fronts recedes to z(s, t) — —zoo (hence Z ~ —i). 
For z and z' along such remote fronts, say located about 
2/ = y < 0, -r > L, the quantities ±iZ, Z/Z' and Z'Z' 
go to ±1 up to exponentially-small terms and all the 
above IC{-,q) functions decay to zero like exp{TrY/2L): 
equation (|5D|) then resumes its free-flame form (|T?)l . It 
also formally does when Im(z)/_L > gets large, since 
Z+l/Z then decays like exp(jz7r/2i), whereby the direct 
influence of the obstacle disappears there (see subsection 
HID. 2). 

D. Numerical flames 

1. 'Kinematic' flash-back 

In section 3 a few steady flame shapes in straight chan- 
nels were determined along with the corresponding effec- 
tive propagation speeds Vp- It is (intuitively-) clear that 
a steady flame cannot settle about the obstacle when the 
injection velocity has Vin < Vp, in which case a flash-back 
of purely kinematic origin would take place, the average 
front location Y{t) along the channel axis eventually re- 
ceding towards y — —oo like {Vin — VF)t < 0, by ([TU|: 
studying this phenomenon might be useful in the context 
of safety. Admittedly, the flash-back of real flames con- 
fined in channels with material walls can be a more com- 
plicated phenomenon than simply based on front kine- 
matics and ideal flows: mixture slowing down in viscous 
boundary- layers and/or interplay with conductive wall 
heat-losses, and chemical heat-release often contribute to 
induce it [39[; yet simpler experimental situations, where 
lateral wall effects are not crucial, often exist [io^l, when 
the incoming fresh gas has a flat velocity profile (like 
here). Next, for fiames kinematically positioned ahead of 
a central obstacle at a distance larger than or comparable 
to C = 0(15-^), one may omit direct counter-streamwise 
conductive/ viscous transfers to/from the obstacle, since 
those decay exponentially quickly over an 0{£) distance 
to the latter. One may also note that the obstacle and 
its images in the channel boundaries constitute a 2L- 
periodic configuration, whereby the x = ±L lines need 



not represent material surfaces, only lines of symmetry, 
and hence do not necessarily bring about viscous effects. 

In the ideal-flow simplified formulation we consider 
here, the way the a:-averaged front location, y < 0, varies 
at steady state when Vin — Vp = SV approaches zero can 
be estimated as follows, if one anticipates that the front 
shape can be described by y = F -I- F{x) -t- 6F{x): here 
F{x) is the steady free pattern belonging to Vp (prop- 
erly shifted, since Ea. ([n)) is translation- invariant along 
the y-axis) and SF{x) is a small correction to its shape 
induced by the obstacles at large, but finite, distance 
from the front. As seen from (j3Gp the obstacle- induced, 
explicitly written, terms vanish as Y — — oo; more pre- 
cisely they vary asymptotically like e = exp{TTY/L) <C 1, 
because the grouping Z -\- 1/Z — 2/ sin(7rz/2L) and 
the coefficient-functions accompanying it as factors are 
0{y/e) for such remote flames (see the remarks below 
pop). In the same limit other terms also appear when ex- 
panding n-u\^—n-Uo\^ : these are 0{SV) , or are linear 
functionals of SF(x) with 0(1) kernels (e.g., coming from 
expanding Wo{z, z'), or the metric factor ds/dx); besides, 
pHl) indicates that 6V = 0{dSF{x)/dx) at steady state, 
since 2LVin and 2LVp both measure front lengths. One 
is led to conclude that SV and 5F{x)/L generically are 
of same magnitude and, if the obstacle is to have an in- 
fluence on y, both need be 0(e). This simple order-of- 
magnitude balance leads to: 

T:Y/L = \n{{Vin-Vp)/a) <0, < a = const. (31) 

for Vin just above Vp. As shown in the next sub- 
section, Vin ^ 1 leads to 0{VinL) flame heights hence 
< y = 0{VinL), suggesting that the x-averaged flame 
location monotonically increases with Vin'- in accordance 
with pip for Vin in near-flash-back conditions, and lin- 
early for V^n ^ 1. 

At least the above is one's expectation when the steady 
free solution is unique. In case of multiplicity of steady 
g = profiles, each with its own Vp, or if the barely set- 
tled fronts happen not to stay steady, the situation is less 
clear. As pointed by an anonymous referee, the existence 
of multiple solutions indeed raises questions of principle 
as to 'the' flash-back limit even in steady noiseless situa- 
tions, e.g., when Vin is intermediate between two allowed 
steady values of Vp, especially if these are nearby. If 
the profile belonging to the largest speed Vp = Vp'^ has 
a special symmetry and the generic one has the second 
largest Vp = Vp^ with Vp^ > V^n > Vp \ a transient 
flash-back may take place until the symmetry is destroyed 
by the slightest disturbance, leading to a front that ulti- 
mately settles near the obstacle. This can become really 
complicated if the very fact that the front approaches the 
obstacle restores the aforementioned symmetry (hence 

Vp ~ Vp^ > Vin), as seems to be the case in Fig. Ufa): 
oscillating Y{t) are not excluded then. In the simpler sit- 
uation implicitly assumed to obtain pip, where max(yp) 
belongs to the generic asymmetric pattern like in Fig. [2j 
the above estimates carry over to unsteady (yet quasi- 
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FIG. 4: a (left) Numerical fronts near an obstacle for E — 5, 
L = 2, £ = 0.1 and, from bottom to top, Vin ~ 1-5, 2.0, 3.0 
and Vin ~ 10 (moderate spatial resolution), b (right) Same 
as for top curve of (a) with a finer resolution, and associ- 
ated streamlines (thin curves) refracting across the front; the 
dashed line is the analytical result, Eq. ()33|) with Sy/L = 0. 



steady) situations, provided Y{t) ~ Y = 5Y = 0{L). 
From (fTU)) they lead to 

^ = M„-^^i=^-aexp(^(y + 5r)) , (32) 

with the same constant a and Y as in pip : as an- 
ticipated from (Vin — VF)St ^ 6Y, the natural time 
scale diverges near Vin = Vf^ t = 0{L/e), and flash- 
back is to be sensitive to 0{e) fluctuations (of what- 
ever origin) SVp^t) = Vpit) — {Vpit)) about the time 
averaged flame speed {Vpit)). In case SVpit) is ran- 
dom and evolves on the t — 0{L) scale, as happens 
with and the ZT equation in the presence of small 
noise, ([5^ would become a Langevin equation for the 
flapping 'motion' of SY{t) in the potential —SYiVia ~ 
(Vf)) + (aiexp(7rr/2i)/7r)exp(7r5y/L). A smaller pos- 
itive Vin — (^f) clearly means a less confining potential in 
the negative SY direction, and would result in larger fluc- 
tuations in mean flame location Y ^SY along the channel 
axis for a given diffusion coefficient associated with the 
variance (((5yF(i))^) of the fluctuations and their time- 
scale. But how large/fast really are such fluctuations? 

To take up the question numerical runs were devoted 
to a few Vin about Vp, producing the three lowest curves 



in Fig.Ha) (E = Z = 0.1, L = 2, g = 0.1). Those 
two in the middle are steady fronts with F(~x) — F(x) 
and correspond to V^in = 2 or Vin = 3 (from bottom to 
top). The lowest front has V^n = 1.5 and is likely headed 
to y = — oo (recall that Vp — 1.58 or 1.47, or... smaller. 
Sec. 3). Its shape is still unsteady, however, the extra 
crest it exhibits being headed to a; = +L, and another 
one is to crop on top of the left trough; nothing indi- 
cated that the widest profile with a single trough like in 
Fig. [5] will ever be recovered. The issue partly depends 
on numerical noise: quasi-steady fronts rapidly become 
sensitive to it as their wavelength (A = 4L > 25A„ in 
Fig.m^a)) increases [s^, and develop sub- wrinkles which 
make Vp, then Y{t), fluctuate. 

Clearly, a fuller study of the flash-back phenomena, so 
sensitive to fluctuations and presenting a critical slowing- 
down, would not fit in here and will be developed else- 
where. To conclude the subsection we nevertheless recall 
that the arguments leading to (|5T|) basically rely 

on a simple order-of-magnitude balance and on three ro- 
bust properties: (i) the infiuence of a 2L-periodic array 
of remote obstacles on the bounded potential fresh gas 
flow decays like exp(27r?//2i) for large negative ys, which 
brings about an 0(exp{TrY (t) / L)) inhomogeneous forcing 
term (~ cos(ttx/L)) in an otherwise free front dynamics; 
(ii) the overall mass balance (jlOp holds; (iii) the fluctua- 
tions in x-averaged front location Y{t) occur over a longer 
time scale than 0{L), which is in fine under control of 
Vin — (Vf)- a simple equation like ((32|) is thus expected 
to hold even if vorticity effects were retained, these essen- 
tially affecting the unspecified coefficient a > therein. 
The problems relating to the multiplicity of nearby meta- 
stable steady states when L 3> A„ will remain, however, 
since they already exist in the Sivashinsky limit E 1^: 
augmenting (jl6p or its ZT generalisation so as to include 
the 0(exp(7rF/i)) residual influence of a remote obstacle 
could constitute a good starting point to take up them, 
because the pole-decomposition method then gives inde- 
pendent access to noise- free steady solutions pl| . 



2. Larger VinS, noise 

As V^n increases beyond 3, the front gets more and 
more markedly V-shaped while staying nearly symmet- 
ric. For Vin S> 1 its height is expected to asymptoti- 
cally scale like LVin since it length does, by ([TU|. Such 
slender steady front shapes could be deduced from the 
integro-differential equations ([T4|((30)) on exploiting two 
simplifications: (i) for < y < 0{LVin), {Z + 1/Z) = 
0(exp(— 7ry/2L)) ^ 1 is neghgible in ([50]) . implying a 
vanishing direct influence of the obstacle which merely 
keeps the flame near x/L = when y/L = 0(1); (ii) a 
4L-periodic linear array of charges soon looks uniform as 
distance to it grows, whereby the kernel Wo{z, z') in the 
ordinary integrals of (jl4p is piecewise-uniform and pure 
imaginary when |y — y'\/L ^ 1 (indeed, tan(a -I- ih) ^ 
i sign (6) for |6| ^ 1). This ultimately reduces the prob- 
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V,^ + (E-l)y/L + 1 



FIG. 5: Tangential velocity u ■ t vs. altitude y, for the same 
flame as in FigUfb). The solid line is from numerics, the 
dashed and dash-dotted lines are u ■ t = Vin + {E — l)y/L + 
0(1), and the markers at \y\ = r = iqL/n locate the obstacle. 

lem to an ordinary differential equation for the flame 
shape F{x) or its inverse x = X(y). 

This can also be obtained directly from (HJ® in an 
easier, more physical, way on exploiting the scale dis- 
parity x/L = 0(1) <C y/L = 0{Vin) over most of the 
region where the flame lies (see Appendix). The analysis 
involves v{x,y) ~ u ■ t = Vin + {E — ^)y/L + 0{l) as 
intermediate formula, and finally results in: 

= i+{E^i){l-\x\/L) + 

where 5y — 0{L) is a constant shift whose exact value 
depends on the structure of the y/L = 0(1) region 
about the plane y = and, in particular, on the ob- 
stacle/flame interactions. Comparisons of the predicted 
(and suitably shifted-) shape, Fig.|4l[b), and tangential 
velocity, Fig.[5l with numerics yield very good agree- 
ments even for moderate values of Vi„, with expected 
magnitudes {0{1/Vin) and 0{1/Vi^)) for the fractional 
errors on u ■ t and its derivative. Since ds/dy — \ ~ 
{dX/dyY/2 = 0{l/VinY < 1: one may identify y and 
the arclength s, with s = at the flame base, where 
(l33t holds. Interestingly enough, (p3)) predicts a leading 
order flame height F{L) — LVin that does not depend 
on the expansion ratio E. The time of transit along the 
front, ^transit, obtained on integrating dt — ds / u ■ t from 
< s 0{L) to Smax = LVin (sec below Eq. ^0^) is 
^transit = L \n{E) / {E - 1) - O {L / Vn) < L; only ioi E ^ I 
(and Vn ^ 1) does this coincide with the value(= L) per- 
taining to unconflned fronts in fast flows. 

The above analysis and result are only outer ones, and 
cease to be applicable in fully two-dimensional regions at 
0{L) distances from the obstacle altitude and from the 
exact flame tip. Fortunately, at such scales, the two sides 
of the front are distant by 0{L/Vin) only and act as uni- 
form line sources of burnt material along the walls (about 



the tips) or X = axis (near the flame base). An analysis 
of such regions could be performed to find the fiow fields 
that match the region 0{L) <y < F{0) - 0{L) with the 
downstream region F(0) < y — 0{VinL), where u = 
and V — EVim or with the upstream one (u = 0, v = Vin)- 
The resulting 0(1) velocity variations would get singular 
near the tips and the obstacle, however, thereby necessi- 
tating still other inner zones there. . . ; One of the virtues 
of integral equations like (ITU) ([501 is that the above nested 
asymptotic structures are all accounted for at once; by 
the same token, they require a definite skill to be han- 
dled analytically without a priori knowledge about the 
location/size(s) of the various zones. 

These asymptotic trends could not be accurately 
checked in steady cases up to very large VinS for two 
reasons. First, as Vin gets high, the front curvature 
at its base, C{x — 0), approaches that of the obstacle, 
1/ry « ^/qL, and the adopted local burning speed 1 — CC 
decreases accordingly; it may even vanish there iiry < C 
and further increasing Vin would result in flame blow-off. 

Before such high velocities were reached, however, nu- 
merical noise came back into play. Using a higher spa- 
tial resolution to reduce it is costly, as the CPU effort 
for the simulation of a flxed time lapse grows with the 
number N of numerical nodes like N'^ {N'^ from the con- 
straint 2C5t < ((5s,nin)^, another factor of N'^ from the 
non-convolution integrals, Eq.([5])), but this suppresses 
the sub- wrinkles; see FigUJ^b), computed with N — 1024. 
The plot also shows the corresponding streamlines ob- 
tained here by integration oi dtz = wiz) (since the front 
is steady) with w{z) given by (15)) ([^ (pSj) : note their re- 
fraction as they cross the front, caused by n • it having a 
jump and t u being continuous there, see ([3]) and above. 

Although the precise origin and way of action of such 
a numerical noise are not under full control hence not 
completely clear, the interpolation jitter caused by dy- 
namical node monitoring where the front curvature gets 
high (near crests or the obstacle), undoubtedly plays a 
role (finer resolutions remove it), and acts as a spurious 
random source of normal node displacements, hence of 
front deformations. The mechanism of amplification of 
small noise- induced wrinkles is the same as in [42, 43j , 
and involves the interplay of curvature effect, hydrody- 
namic instability, wavelength stretching and, here, in- 
teractions with front images. When the tangential ve- 
locity It • T is high and uniform along the unperturbed 
front the instability may only be convective (as op- 
posed to absolute) [l^], and the drifting sub- wrinkles 
triggered by tiny numerical noise would eventually have 
max(ti7(A:))=(£'— 1)'^/16>C as Lagrangian growth rate and 
2A„ = %ttC/{E - 1) as final wavelength (= 0.628 for 
i? = 5,£ = 0.1)if raJmax ^transit is large enough, which is 
experimentally known [42l. |43|. Here, confinement makes 
U T increase along the y-axis (Fig. [5]) and the wavelengths 
A also increase as the sub- wrinkles approach the flame tip 
[iaj . More precisely, since tt-r Vin + {E-l)s/L, Fig.lHl 
the wavelength ultimately increases with s according to 
the "node conservation" law X/u ■ t = const, and hence 
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also ends up varying linearly with current arclength s; 
this can be detected in Fig.HJ^a). By the same token, the 
wrinkles also get closer to the lateral walls, which tends 
to quench their growth in amplitude (not in wavelength) 
through potential interactions with their closest image 
in the walls, even in absence of significant nonlinear ef- 
fects: the \ow-k growth rate w{\k\ ^ /c„) of disturbances 
about a front nearly parallel to an impermeable wall in 
the fresh gases, and lying at a distance A > to it, is in- 
deed reduced by an extra factor of (1 — exp(— 2|fc| A)) [45l |. 
This follows from the fact that a near harmonic wrinkle 
of local shape A{s) exp(z k{s') ds), with the amplitude 
A{s) and k{s) varying over some scale A ^ l/|fc|, induces 
in the fresh gas potential flow a disturbance in u. • n sa u 
that locally decays with coordinate S normal to the front 
like sign(,5)|/c(s)| A(s) exp{~\k{s) S\ + i f k{s') ds')- the 
image 'beyond' the wall at a distance 2A(s) does the 
same... up to a reversed sign of A(s) and the change of 
6 into 2A(s) — (5, which ensures that the slip condition 
u = at the walls is satisfied and does reduce the it- 
fluctuation felt by the front at S — 0. The growth rate 
vu{k) is replaced by: 



7y{k) = 



E- 1 



1 



|fc| , 
— - exp(- 



-2|fc| A) 



(34) 



< 1, 



to be used later (Section V). As soon as 2Ak„ 
Wy{k) < whatever |fc| > 0, implying damping. 

Using the law of wavelength variation X/u- t — const, 
mentioned earlier, one can extrapolate the observed 
A(?/tip) ~ 0.75 (~ 2.4 A„ for the parameters used in 
Fig.m^a)) back to the point y* where u ■ t reaches a local 
maximum, m-t « 16, Fig.[5l giving a \{y*) of 0.25 < A„; 
this is not incompatible with a numerical noise originat- 
ing from a region of size comparable to the obstacle radius 
r = AqL/n, r w 0.255 here. Note that the deformations 
of the left and the right front wings are little correlated 
(if at all), which comforts their interpretation as noise- 
induced. 



V. FLAME SETTLED NEAR LATERAL BUMPS 

Whereas the centred obstacle clearly relates to flame 
stabilisation by bluff bodies in high-speed flows, there 
also exists a 'dual' configuration of interest (theoreti- 
cal, and relating to combustors), where flame stabilisa- 
tion is achieved by means of 'bumps' protruding inward 
from otherwise-straight lateral walls. There exist several 
mathematical ways of 'growing' bumps out of smooth 
surfaces [i^, but here we chose to exploit the material 
already presented above: the selected bumps will consist 
of the two mirror-symmetric halves (left and right) of the 
preceding obstacle, stuck on the walls with the convex 
parts facing inward (Fig. 15]). The first step towards the 
corresponding Green's function again is to conformally 
map the physical plane, now called z, to another one 
where G is accessible. To this end we invoke the supple- 
mentary complex potential ^'(Z, g) defined by (I^Hl) . and 
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FIG. 6: Fronts positioned near lateral bumps, for Vin = 5 
(lower curve) or Vin = 10, both with L — 2, E = 5, C — 0.1, 
q = 0.1. 



notice that the streamlines it generates are symmetric by 
X — ^ — X and x — i— 5-L — x, is regular at finite distance 
along the lines x — ±L, and can be analytically continued 
periodically along the x-axis; this, actually can also be 
seen from the identities (P5])(j26D . So, once properly nor- 
malised and translated, the above 'inert fiow' '^{Z{z),q) 
can provide one with a natural change of complex plane. 
The Boussinesq-like mapping [47| 



c 

dC 

dz 



L 



(35) 



2iL P{iZ{z),q) 
~ P{-tZ{z),q) 

l_l[/C(-zZ,g)-/C(zZ,g)](Z+|) (36) 



between the physical z plane and a straight channel in 
the auxiliary C-plane does the job, whereby the required 
Green's function in the present case merely is: 



G(i,z') = Go(C(^),C(^')), 



(37) 



where G'o(-, •) is the q = Q function defined in (|TT|l . 

The supplementary flow needed in W has in fact al- 
ready been found and, by construction, is simply given 
by Wsup = d^sup/dz, with 



-*C(z) 



E- 1 

4:L 



<j{s') ds' + Vh 



(38) 



IS 



once ([9]) is accounted for. The contribution Wsot 
accessible from ([51) (|37l) . and Wsup on z-differentiation of 



13 



, which ultimately yields: 



n - =Re^—i n{s) ^ /• ( T4n+ / ^ o'(s) ds 



E-1, 
'~8L 



-Re 



g l'n{s)Wo{C,CMs')ds'-ALais) 



(39) 



thereby giving one access to an evolution equation again 
formally identical to ([T5|) 

Whenever the front shape are steady and symmetric 
by i — )■ —X, they coincide with those obtained in Sec. 4, 
up to an horizontal shift by ±L and a mirror symmetry; 
these are not duplicated here. 

The response to small-scale random noise is different. 
In Sec. 4, the sub-wrinkles observed near the flame tips 
were constrained by the nearby walls to form 'varicose' 
type of patterns |45| . i.e., symmetric under a: ± L — >■ 
L ^ X. When the front leading edges sit near lateral 
bumps, however, the flame tip is not that constrained any 
longer and can preferentially develop. Fig. El the 'sinuous' 
(anti-symmetric in x at 0(A„) scale) mode of wrinkling 
because it has a larger growth rate than the 'varicose' one 
(ill |45| : the 'DL' effects of the near-facing sides of the 
front, at a mutual distance of 2A > 0, then reinforce one 
another instead of mutually cancelling partly as they did 
for the varicose mode (see end of Sec. 4), which enhances 
the low-fc growth rate of harmonic disturbances by an 
extra factor of (H-exp(— 2|fc| A)) [ill] whereby the (again 
Lagrangian) growth/decay rate zu becomes: 



TUs{k) 



E-1 



\k\il 



exp(-2|fc| A) 



(40) 



to be compared with ^^{k) in (1341) . Note that vos^max, 
reached at k = k„ and A = 0, is 4cLii,^,„ax = 4uj„iax, im- 
plying a strong preferential amplification of the sinuous 
mode near the flame tip. 

Both growth rates can be used to estimate the ra- 
tio in final sub-wrinkle amplitudes with interaction be- 
tween near parallel fronts {Ay or As). Instead of lineariz- 
ing (|39p about a symmetric steady pattern y = F(x) - 
then trying to solve the linear problem at the expense 
of extra heav y a lgebra, a simple WKB type of physi- 
cal argument [4^ is presented below when T^n ^ 1 and 
k ^ km kL ^ 1. To this end we again interpret Ws and 
Wv as describing disturbance growths with Lagrangian 
time T(s), related to s by dT = ds/u ■ r [i^. Assum- 
ing — 1 when the unperturbed front sides are 
far apart, the final ln|As/A^| is obtained on integrating 
{ws-vjy)dT = {E~l)\k{s)\ exp(-2|fc(s) A(s)|)ds/M-T 
along the entire front. The needed u ■ t and A(s) = 
|X(j/)| follow from the analytic steady shape y — F{x) 
at T^n 1 (or its inverse x = X{y)), as is deduced 
from the properly shifted Eq. ([55)1 (|x| replaced by L — |x| 
there) and the corresponding linear variation of it • r; 
fc(s) = 2it/X{s) is given by the 'node conservation' law 
X{s)/u ■ T — const. Acknowledging that the dominant 
growth of As/ Ay takes place near the flame tip, one may 



even approximate {zus -— Wy) and A by their late be- 
haviour (i.e., set A = |s — StiplZ-EV^n and freeze k{s) to 
some fco ~ kn therein). The final (fco-independent) esti- 
mate of the sinuous-to- varicose amplitude near the flame 
tip reads: 



=exp((£;-l)/2). 



(41) 



With E = 5, this is only 7.34 but sufficient to make 
Ay still hardly visible when the sinuous mode becomes 
detectable {Ag = A/10, say) near the flame tip. Fur- 
thermore, the newly born cells are too curved once fully 
formed to be further affected by noise [H, H^] before they 
reach the tip, and in no way can the varicose mode emerge 
once the sinuous mode has appeared. Put in words: the 
wrinkled front responds to random noise as to avoid self- 
intersections hence the formation of blobs of fresh gas 
detaching from the main front's tip; this is experimen- 
tally known [i^] and intervenes in how to reduce emis- 
sions of combustion- induced noise [i^. Like in Fig. |3] the 
sub- wrinkles noticeably increase the front length per unit 
length as they travel up the y-axis, which gradually re- 
duces the total flame height since the time-averaged front 
length is 2L14n in any case, by p^ : this flame-brush 
shortening by cusped cells is well known experimentally, 
even in some turbulent flames [13, • 

Very large amplification ratios proportional to 
exp(ti7niax^transit) explain why even residual, hardly mea- 
surable, turbulence is enough to trigger the appearance 
of sub-wrinkles along unconfined Bunsen flames at high 
ambient pressure po [i^: simply because the grouping 
■ci^maxitransit scales like y/po for most usual hydrocar- 
bon/air pre-mixtures (C/q ~ I/Vpoi -^th ~ 1/Po), which 
results in huge variations of its exponential as po increases 
(up to 100 bars in [1^), at fixed VSn/C^o not to modify 
the flame height-to-width ratio. Possibly for related rea- 
sons, many laboratory burners tend to have compara- 
tively small values of vumux ^transit , resulting from rather 
small L/Xn ratios, smaller than the L/Xn ~ 6.5 in Fig.^ 
the little confined Bunsen-burner flame in figure 3 of [50| 
has L/Xn ~ [Xn ~ 100 £ for lean propane/air mixture 
like the upper half of the Vm = 10 front shown in 

Fig.m 



VI. BREATHING FLAMES 

An evolution equation for the flame front (only), with 
built-in boundary conditions and each operating param- 
eter adjustable at will, constitutes a flexible tool. The 
example considered next to make the point, selected be- 
cause it relates to flame response to incident pressure 
waves [13, [m , deals with fronts near an obstacle then at 
lateral bumps, and fed at an oscillating injection velocity 

Vi^{t) = Vin, av + 5Vin sin(wi„t). 

Before proceeding to the results of numerical integra- 
tions proper, it is of interest to again consider the large- 
V^n limit, now in unsteady situations. The analysis lead- 
ing to (j33p can indeed be adapted almost mutatis mu- 



tandis to slowly- varying flames and velocity components; 
this merely requires to include V — dtX/{l + X^)^/^ 
dtX in the kinematic conditions (|A1|) and gives a partial 
differential equation for X{t^s): 

dtX + ds{X[{E - \)s/L + = -1, (42) 

in which the inlet velocity Vln is assumed to evolve on the 
t = O(ttransit) time scalc. The hyperbolic Eq. (|42|) can 
be solved with the boundary conditions X{t,6y) = o{L), 
corresponding to a front leading edge near an obstacle 
for Vin ^ 1, on the grounds that the left undetermined 
Sy = 0{L) shift in (|33| does not depend on T4n, and hence 
is time independent here; a i-dependent 6y would shift 
V;„(i) to T4n(i) - d5y/dt in (gS))- The long-time solution 
to (HI]), once written in the form y = F{x,t) like in 
reads 

F{x,t)-Sy^ / exp{{E-l)t'/L)Vinit-t')dt', (43) 
Jo 

with (£; - l)e(x) = -iln(l - {E~ l)\x\/LE), so that 
6(0) = and Q{L) — itransit- For Vin = const, this 
resumes and when Vin(i) is periodic, the time- 

averaged {F{x,t)) is still given by ([55]) with Vin replaced 
by {Vin{t)). With the aforementioned oscillating Vin{t), 
(|43p produces wavy fronts undulating in a harmonic way 
and lagging behind the injection velocity with a phase 
shift encoded in Q{x); these are not unlike either side of 
(nearly-) confined experimental Bunsen flames fed with 
a slowly osciUating Vin (e.g.. Fig 9(b)(c) in (s^). Yet 
there exists an important difference caused by the more 
elongated shape, relating to the pulsation needed 
to "imprint" A^ceii > 1 cells along a fiame of half front 
length LVin- For unconfined or weakly confined Bunsen- 
or V-flames [H, \5^, one needs Win = 2T:Nce\\/L (= 
15.7 for L = 2, A^ceii = 5); with the present confine- 
ment, our model requires a pulsation larger by the factor 
L/ttransit = " l)/ln£; > 1 (= 2.5 for £: = 5); a close 
scrutiny of the E = 1 and E ~ 5 curves in Fig.[71[a) makes 
the point. Once combined with the related change in 
steady profiles, also caused by the confinement, this vari- 
ation in cell number and the accompanying wavelength 
stretching quite logically are somewhat at variance with 
experiments on the usually unconfined (4^ , or barely con- 
fined [i^l, burner flames. Another noticeable feature of 
Eq. ([1^ is the absence of any DL-like mechanism of in- 
stability at this order: a damping term —{E—l)X/Leven 
appears on expanding the braces of Eq. (j42p . stemming 
from a negative (dxu)- ^ ~dyV < and thus resulting 
from the confinement that breaks the translation invari- 
ance along the x-axis. To explain this absence of any DL 
mechanism in Eq. (j42l) we note that the wave number 
belonging to such wavy fronts is 0(l/LVin), whereby the 
associated growth rate w{k) ~ |A:| is negligible compared 
to 9tln|X| = 0{1/L). Such a conclusion would thus 
still hold true even if the flame were not symmetric any 
longer... yet still varying over the same long time-scale 
as above. 
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FIG. 7: Breathing flames in the presence of (a) An obstacle, 
with Vin.av = 10 and E — 5, SVin = 6 (thick line) or i? = 1, 
SVin = 10 (thin line); (b) Bumps for E ^ 5, SV,n = 4.5, and 
Vin,av = 6 (lower curve) or 10 (upper curve). In all cases 

CC7in = 40. 



The hydrodynamic instability will thus manifest itself 
only at high enough pulsations Win that the resulting cells 
along the front have A < 0{L), which here requires higher 
WinS: to fix the idea, the preceding estimate with the in- 
dependent factor of 2.5 included, gives Win — 39.3 for 
A'ccii = 5. In turn this will bring back the nonlinear and 
curvature effects (also negligible in (|12|)), thereby neces- 
sitating numerical integrations. Note that the natural 
exponential decay of the u-disturbances with distance to 
the flame and its images in the walls (see above equa- 
tion ([51)1 ) strongly reduces the flame/images interactions 
when A < 0{L), thereby somehow restoring translational 
invariance for the shorter wrinkles (especially when ap- 
pearing near the tip and along the channel centreline); 
in particular, these get little subjected to the damping 
featured in (|42)) . 

Figure[7]shows numerically obtained fronts settled near 
an obstacle or bumps, fed at the same oscillating veloc- 
ity Vin{t) = T4n,av + SVin sm{ujint) , with L = 2, C = 0.1, 
q = 0.1, Win = 40, and E = 5 in most cases (see cap- 
tion of Fig.[7|). To guide the eye, a front corresponding 
to the same parameters as above but without gas ex- 
pansion {E — 1 in lieu of = 5), is also plotted as a 
thin line in Fig.[71[a): suppressing the DL-like instabil- 
ity mechanism for frequencies leading to A < 0{L) does 
affect the amplitude of the flame response, despite the 
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larger SVin = 10 instead of 6, and gives a larger Nccw 
because this also increases the itransit/i ratio to 1; this 
also suppresses the wavelength stretching. As was the 
case for the growth rates w{\k\ ^ fc„) then about the 
effective speed Vp of free flames, the present formula- 
tion plausibly overestimates the consequence oi E ^ 1 
on the amplitude of wrinkling, however, especially near 
the tip(s). This trend of potential-flow models should 
be kept in mind when comparing the present results at 
E 1 with data on actual flames responding to inci- 
dent pressure waves; yet flames anchored at the rim of 
a Bunsen burner look like Fig.[71[b) when excited by a 
loudspeaker to generate oscillating VJnS [5l|... once the 
selected frequency Win properly accounts for the shorter 
time of transit along the unperturbed front when a con- 
finement is present. 

One can note the appearance of a sinuous mode of 
sub- wrinkling at high mean injection velocities, Fig.[71[b), 
again triggered by numerical noise: if allowed by the 
problem symmetries, this eventually catches over the 
varicose-type of mode generated by a pulsating uniform 
injection, and leads to flame tip flickering. Although not 
unphysical (unwanted flickering often appears in experi- 
ments) and removable here by using a finer - yet 'costly' - 
resolution, this influence of noise should be under better 
control and accounted for explicitly, instead of appearing 
as a side effect. The implicit method of [52,] shows good 
prospects of improvements because it allows for much 
larger time-steps without sacrificing spatial resolution. 
However, this will not change the fact that each of the N 
nodes in the discretized integral evolution equation(s) is 
coupled at any one time step to its TV — 1 neighbours, in 
the long-range manner that characterises d = 2 Laplace 
problems, whereby evaluating it • n|_ at time t requires 
0{N'^) operations; adapting hierarchical algorithms de- 
veloped for similar A^-body problems (e.g., [2^, [s^) to 
reduce this to O(AlnA^), yet at the expense of high 
noise levels, would greatly help fully exploit the evolution 
equation approach, at least for flames randomly forced 
on purpose where statistics (hence long runs) are often 
needed. 



VII. CONCLUSION, PERSPECTIVES, OPEN 
PROBLEMS 

In three conflgurations (free evolutions, obstacles of 
two different kinds) we showed that Green's functions 
and conformal map methods, combined with numeri- 
cal integrations of the corresponding Frankel type evolu- 
tion equations we derived, can indeed simultaneously ac- 
count for flame/flow nonlocal interactions, large front de- 
formations and nontrivial boundary geometries (includ- 
ing multiply-connected domains). Once the numerical 
method is properly optimized, this approach will be able 
to handle flames of even larger lateral extents 2L than 
here (in Fig.|4j 2L 1300^ « 8 cm, for the same mix- 
ture as in figure 3 of [soj). 



The shapes/evolutions of channelled fronts so obtained 
share a lot with actual flames: hydrodynamic instabil- 
ity, formation of cusped cells, sensitivity to noise when 
wide enough (or at high pressure), sinuous mode of fi- 
nal wrinkling, asymmetric free flames; yet the potential- 
flow models were found to overestimate the consequences 
of density changes on the wrinkle rate-of-growth and fi- 
nal amplitudes, while the laboratory fronts seldom are as 
confined as the channelled ones considered here - with 
identified consequences on the flame morphology - and 
often are of smaller lateral extent. 

Besides the above mentioned physical and numerical 
issues, the present models raise interesting mathematical 
related questions, two of which are now evoked. First, 
it would be interesting to explore whether potential-flow 
models of steady flames with E — 1 — 0(1) admit solu- 
tions in the limit v = nC/ (i? — 1)L ^ 0+. And if they do, 
like (ITSl) (1171) , is there again a continuum of allowed flame 
fronts 34] that is broken into a densely-packed discrete 
set thereof (with very nearby Vp\) when small-i^ effects 
are restored (which relates to flash back)? Second, as the 
steady version of ^ with a prescribed potential fresh-gas 
flow- field n • m|_ — V(/3_ and the Poisson equation with 
given sources along a prescribed front y = F{x)) [sol. 
admit separate variational formulations, is there an ex- 
tremum principle for steady potential-flow flame mod- 
els? (And, if 'yes', what does the corresponding 'action' 
physically mean?). This could possibly help clarifying 
why potential-flow apparently overestimate the influence 
of i? — 1 > in comparison with actual flames, in case 
the 'action' in question would have a second variation of 
definite sign. 

One must indeed not forget that actual flames are also 
affected by vorticity, il. As regards it, flrst notice that 
the expression for the complex velocity w at the entrance 
of the free flame fronts studied in Sec. 3, Eqs. (O (|12p . 
can be re-written in operator form as: 

2i ^"^"^^ 

where [n • it] = {E — 1) a{s) measures the strength of 
volume sources along the front, the operator j7o is minus 
the principal-part integration over ds' Wo{s, s')/4:L ap- 
pearing in (1121) , and I is the identity. The last equation, 
(|44t . closely resembles Eq. (42) obtained in jl^ for steady 
flames in straight channels, and is actually an integrated 
vorticity-free version thereof. The integration constant 
is determined in the present = model via an over- 
all mass balance; on the other hand, the ^-differentiated 
equation in |18i] has a distribution of vortices along the 
front besides (modified) sources of volume, which is miss- 
ing here and formally amounts to supplementing ([3]) with 
a suitable $7+- and ?«_-dependent jump [m • t]. As shown 
above, the shape of non-straight channels or obstacles 
can be encoded in the Green's functions, which here 
changed the Jo operator into a more general form J7g, 
since G{Z,Z') ^ Go{Z,Z'), see dST]); this also made the 
right hand side of (|44l) a known function of z(s) stemming 
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from Wsup ^ Wsupl— «oo), i.e., from the nonlinear maps 
from the physical to the auxiliary planes (for example, 
see Eq. (l38l) '). which can thus be adapted to more general 
channels than the asymptotically-straight ones consid- 
ered here. The important step to take up now clearly is 
to investigate how the approach of [l^ adapts to such sit- 
uations. After all, the Green's function associated with a 
point-vortex in the annulus of Fig. [3] and Neumann con- 
ditions also is expressible in terms of the same function 
P{-,q) as in (1^^ . This will be pursued elsewhere. 

Appendix A: Slender steady flames (central 
obstacle) 

Thanks to the scale disparity x = 0{L) <^ y — 
0{VinL), dxx ^ dyy in the Laplace equations sepa- 
rately satisfied by u{x,y) = 0{1) and v{x,y) = 0{Vin), 
whereby both functions are affine functions of x with 
y-dependent coefficients, to the two leading orders in 
l/14n ^ 1- The slip condition u{L, y) = gives u{x, y) = 
a-{y){x-L) = 0(1), v{x,y) = A_{y)x + B_{y) for 
< X{y) < x < L and < y < F(0); u{x, y) = a+{y) x 
and v{x,y) = A+{y)x + B+{y) for < a; < X{y). The 
zero-vorticity condition dxV — dyU — 0(1/Mn) implies 
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